Relaxation to equilibrium driven via indirect control in Markovian dynamic^] 
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INTRODUCTION 

In the past decades quantum mechanical systems have 
attracted lot of attention for their peculiar properties, 
that indicate they are good candidates for the implemen- 
tation of outperforming technologies in the fields of infor- 
mation and computation^]] . In this spirit, some amazing 
protocols have been recently developed, as for example 
the computational algorithm for the factorization of a 
large number 0] , or the schemes for teleportation 0] and 
quantum cryptography Q. Many concrete physical sys- 
tems have been proposed for the practical implementa- 
tion of these ideas, as optical devices, cold trapped atoms, 
nuclear spins in magnetic fields (NMR), or quantum dots 
in electromagnetic cavities. 

In all these cases, the largest obstacle to the implemen- 
tation of stable and efficient schemes is represented by 
the unavoidable interaction of microscopic systems with 
the surrounding environment. Because of this interac- 
tion, the system dynamics is subject to loss of coherence, 
irreversibility and dissipation, and the appealing proper- 
ties of quantum systems are usually lost or compromised 
during the time evolution. 

Usually, the environmental action is accounted for 
by describing the dynamics of the system S through a 
Markovian one-parameter family of maps {74 ;t 0}, 
satisfying the semigroup property 7 t+s = Is, with 
t, s 0, with 



Ps(t) =7t[Asr(0)], 



(1) 



where the statistical operator (or density matrix) ps is 
an Hermitian, positive, unit trace operator, acting on the 
Hilbert space associated to the system, and representing 
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its state. This representation of the dynamics is not the 
most general, but it is well justified in many cases, in 
particular when the coupling between S and the environ- 
ment can be considered weak. The generator L of the 
dynamics can be obtained by writing (JlJ in differential 
form, p s = L\ps\, and it has the standard structure 



L[ PS ] = -i[H 8 , p s ] + ]T dj (F iP sFj - l{F}F u ps} 



hi 



cussions. 
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(2) 

where Hg = Hg is the system Hamiltonian, and the set 
{Fi;i} satisfies TrFj = 0, Tr(F t F]) = %. The Kos- 
sakowski matrix C = [cy] must satisfy C' = C ^ in 
order to guarantee the complete positivity of the evolu- 
tion, and then the physical consistency of the formalism 
[E HI ■ It encodes the microscopical details of the inter- 
action between system and environment. The first term 
in the right hand side of (|2|) represents the coherent part 
of the evolution, and the generator of the system dy- 
namics has this form whenever the interaction with the 
surrounding environment can be neglected. The corre- 
sponding time evolution is given by a group of reversible, 
unitary transformations. The second term is responsi- 
ble for irreversibility and dissipation, since it produces a 
contraction map on the set of states, and, in some cases, 
relaxation to stationary states. 

It is of fundamental relevance to study methods to 
fight against decoherence. When the environmental noise 
exhibits some particular symmetry properties, this task 
can be realized by encoding the relevant information in 
suitable Decoherence Free Subspaces or Subsystems unaf- 
fected by decoherence (for a review, see Q). An active 
approach consist in directly affecting the system dynam- 
ics, in order to preserve its relevant properties or induce 
arbitrary manipulations. This controlled evolution is re- 
alized through some functions, entering the dynamics, 
that can be manipulated via external actions (see for ex- 
ample [1, H, [l(| El, E2> EH E3 for a geometric approach 
to controllability). 

Several approaches to the control of a quantum system 
have been proposed in the past years. In the open loop 
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schemes the control functions are a priori fixed (that is, 
they are independent on the state of the system). Con- 
versely, in the closed loop control schemes, the control 
functions are updated in real time by feeding back some 
information about the actual state of the system, usually 
gained via an indirect continuous measurement (quantum 
feedback [IS EE B3] ) • 

The control functions usually affect the Hamiltonian of 
the system Hs, since the environmental action is usually 
uncontrollable. This approach is called coherent control, 
as it affects the coherent part of the dynamics. Moti- 
vated by different experimental scenarios, another con- 
trol scheme has been introduced, in which an auxiliary 
system is used to manipulate the target system through 
their mutual interaction. This indirect control scheme 
is of relevance whenever the system dynamics cannot be 
directly accessed [HI 0> US]. It represents a comple- 
mentary approach to controllability, with interesting fea- 
tures concerning the purification of mixed state s 12 ll , and 
when applied to the dynamics of open systems [22] , since 
it makes use of the correlations between the two subsys- 
tems, that can be created by the environmental action 
(described for the first time in (23jl. 

One of the unwished consequences of the environmen- 
tal action on the system dynamics is, in many cases, the 
collapse of the system into a -in many cases, unique- equi- 
librium state, with a consequent reduction of the reach- 
able sets, and loss of control. In this work we address the 
following question: is it possible to modify the station- 
ary states of a target system T, evolving under a quantum 
dynamical semigroup, by means of an open-loop indirect 
control? In other words, we introduce an auxiliary sys- 
tem, a quantum probe P, couple it to T and consider 
the evolution of the joint system S = T + P, and fi- 
nally discard P by taking into account only the degrees 
of freedom of T. Assuming that S is still described by 
a quantum dynamical semigroup, we study the station- 
ary states of the system T alone, affected by P through 
the correlations between the two systems. The impact of 
both initial correlations, and correlations created during 
the joint evolution, is taken into account. 

The plan of this work is the following. In Section Q] we 
review some algebraic tool for the determination of the 
stationary states of a quantum dynamical semigroup. In 
Section|n]we specify the dynamical settings considered in 
this paper and we derive the relevant algebraic quantities, 
introduced in Section [IJ In Section IHT1 we describe all the 
possible scenarios for the stationary states, in terms of 
the dynamical parameters characterizing the semigroup. 
In Section IIVI we summarize our results, describe their 
physical significance and finally conclude. 



I. STATIONARY STATES OF QUANTUM 
DYNAMICAL SEMIGROUPS 

In general, the second contribution in the right hand 
side of ([!]) leads to the appearance of attractors in the 



state space of S, and consequently relaxation to equilib- 
rium of the states of the system, absent if there is not 
interaction with the environment. A stationary state for 
the dynamics, pf, is determined by the condition on the 
generator = 0. This is a system of linear equa- 

tions that can be solved using standard algebraic tools. 
In this spirit, quantum dynamical semigroups have been 
classified in terms of their relaxing properties [HI]. In 
the uniquely relaxing semigroups, there is a unique sta- 
tionary state, and every initial state eventually collapses 
to it. In the relaxing semigroups, although every tra- 
jectory collapses to a fixed state, this state is not unique 
but depends on the initial conditions. Finally, in the non- 
relaxing semigroups, oscillatory solutions survive. Even if 
this method is very general, the resulting algebraic equa- 
tions are complicate, therefore we will rely on a different 
approach. 

For the Markovian dynamics @ , necessary conditions 
for the existence of stationary states and for the conver- 
gence of ps(t) to them have been derived in terms of the 
operators appearing in the diagonal form of ((2J), 

L[p s ] = -i[H 3 ,Ps] + Y.{ V ^ ~ l{V?Vi,ps}). (3) 

i 

The following theorem summarizes these conditions [1^] , 
and it will be the basis of our analysis. 

Theorem 1 Given the quantum dynamical semigroup 
Q), assume that it admits a stationary state po of maxi- 
mal rank. Defining M. = {Hs, Vi, V}\ %)' , the commutant 
of the Hamiltonian plus the dissipative generators, and I 
the identity operator, the following conditions hold true: 

1. If M. — span(I), then po is the unique station- 
ary state. Moreover, if {Vf,i} is a self-adjoint set with 
{Vi]i} — span(I), then for every initial condition ps(0) 

lim p s (t) = p . 

t — >+oo 

2. If M ^ span(I), then there exist a complete family 
{P n ;n} of pairwise orthogonal projectors such that Z = 
MC\M = {P n ;n} . If{Vi]i} = M., two extreme cases 
together with their linear superpositions may occur. If 
Z = M, then for every initial condition ps (0) 

lim Ps(t) = Y,M P nPs(0)P n ) T ^ p P ° Pn pV 
t^+oo Tr{P n p P n ) 

If Z = M. , then for every ps(0) 

lim p s (t) = VP n /) S (0)P n . 

t— >+oo * — • 
n 

Therefore, in order to characterize the stationary states 
of a quantum dynamical semigroup, it is necessary to find 
a maximal rank stationary state po, and to evaluate the 
algebras M. and Ml '. 
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II. DYNAMICAL SETTINGS AND RELEVANT 
ALGEBRAS 

We assume that S = T+P is a bipartite system, where 
T and P are two copies of the same two-level system, sep- 
arately interacting with a common environment accord- 
ing to the Markovian dynamics (|2|) . The operators Fi are 
given by Fi — Oi ® I for i — 1,2,3 and Fi = I £g> £7j_ 3 for 
i = 4, 5, 6, where I is the 2-dimensional identity operator, 
and <7j, i = 1,2,3 are the Pauli operators. We consider 
the standard representation of these operators in which 
03 diagonal. The matrix C has the form 



C = 



A B 
St A 



(4) 



where A = At is the Kossakowski matrix for the system 
T (or P) alone, and B represents the dissipative coupling 
between the two parties. The form (0J is not the most 
general joint Kossakowski matrix. More complicate ex- 
pressions should be taken into account when three-body 
contributions are relevant (common interactions between 
the two subsystems and the surrounding), and when the 
dissipative couplings of T and P to the environment are 
different (for example, in a non-homogenous medium). 
We will limit our attention to models satisfying (H}; more- 
over, for simplicity, we will further assume B = B< . This 
assumption highly simplifies the mathematical formal- 
ism. 

Following Theorem [TJ we need to write C in diagonal 
form in order to find the operators Vi appearing in ([3]). 
This is achieved by means of the unitary transformation 



UCU't =diag(Ai,i= 1,...,6), 
where Xi are the eigenvalues of C, U is of the form 



y/2 



U U 

-u u 



(•5) 



(6) 



and U, U are unitary transformations such that 
U(A + B)W = diag(A+,i = 1,2,3), 
U(A-B)& = diag(\r,i = 1,2,3). 



(7) 



The eigenvalues of C are ordered as A, = for i = 1,2,3 
and Xi = X~_ 3 for i = 4, 5, 6. Comparing the generator 
forms ([2]) and ([3]), and using the notation U = [iiij], we 
have 



i = 1, . 



k=l 



Following ([6]), it is possible to write 



I ® crj + di <g> I, 



--Vi = 



where we have defined 

3 

k=l 



1, 2, 3 
4, 5,6 



fe=l 



(8) 



(9) 



(10) 



and we used the notation U = [uij], U = [?%■]. The op- 
erators in (fTU|) satisfy Tr<7i = Tr&i — and Tr(<7jcrj) = 
Tr(CTiO't) = Sij. They are self-adjoint if and only if the 
unitary operators U and U are orthogonal. 
The commutant of Theorem [T] can be expressed as 

{H s ,Vi,V?;i\\i^0}' = f| {Vi,V?}'n{H s }\ (11) 

where only non-vanishing eigenvalues Xi have to be con- 
sidered, otherwise the corresponding Vi do not appear in 
the generator ([3]). Moreover, for a given i, 



{Vi,V?}' = {v\ve{ViY,vi sty}'}, 



(12) 



therefore we can limit our attention to the sets {Vi} . We 
find convenient to consider separately the two kinds of 
contributions defined in To begin with, we consider 
a fixed index i such that A^~ 7^ 0, and assume that the 
corresponding di is non-singular. In this case it can be 
written as 



&i — jliRiay,Ri 1 



where 



Ri — Ri 1 — — 



u* 3 + Mi "*i ~ i u i2 
u* iX + i u* 2 -u* 3 - jii 



(13) 



(14) 



and 



Mi 



£(u«) 2 > * = («*3 + Ai)- (15) 



Since I ® ct^ + cr, ® I = jUi72.,(E cr 3 + 0-3 ® I)7£j, with 
IZi = Ri <g) Ri, it follows that 

{I (g> CT, + CTi <g> I}' = <g> (T 3 + cr 3 <g) I}' 7^ (16) 

and then, after the explicit computation, 

{V}' =span(I®I,I®a i! ai<g)]I,CTi(g)CT i ,n + ,A7), (17) 

having defined the additional operators 

fl + = <J\ ® (7i + 02 ® cr 2 + 03 ® 03 , 

Ar = ^(cri®^ -^gxri)^. (18) 

Notice that, in general, the operators in the right hand 
side of (TTT|) are not self-adjoint, nor orthogonal each other 
in the Hilbert-Schmidt metric, since the transformation 
IZi is not unitary. However, if the coefficients u*j, j — 
1, 2, 3, are real, di is self-adjoint and IZi unitary (and self- 
adjoint). Consequently, in this case the basis of {Vi} is 
made of Hermitian, orthogonal operators. 
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Cases 


Conditions 


M basis 


A4' basis 


Z = Mn M' 


{Pn;n} 


I 


n+ = 1, A = A T , B = A 


crj ® CTj, fl+, AT" 


I (8 I, <7j gs 5-j, 
I 8 &i + &i 8 I 


M' 


Pi = TZiHiIZi, P 2 = KiHiTli, 

p :i = 7e»(n 2 + n 3 )% 


II 


n+ = 1, A / A T , B = A 


1(81,0+ 




M 


Pi = ni + n 4 + n+, 
p 2 = n_ 


III 


n+ = n- = 1, A = A T , 
B = aA, a G K n {-1, 1} 


1 8 I, a, 8 5-j, 
I (8 5-j , a; ® I 




M 


Pi = TZiHilli, P 2 = IZiHiTli, 
P 3 = Tl^Ui, Pa = TZiHiTZi 


IV 


n+ > 1, B = A 


I® 1,0+ 


X' D M 


M 


Pi =ni + n 4 + n+, 
Pi = n_ 



TABLE I: Relevant algebras for the determination of the stationary states for the two qubits system, under the assumption 
Hs = 0, in the non-trivial cases. 



The commutants {Vi} are completely characterized 
for i = 1,2,3. Finally, {V,V^}' can be found by con- 
sidering (fT2| : 



{Vi}' 



iff Gi = erf; 



span(I ® I, f2+), otherwise. 

The corresponding sets for i = 4,5,6 can be found 
by applying the same procedure to <7j, assuming that 
\r ^ 0. The result is 

{Vi}' = span(I® 1,1® ct;,ct;® I, a t ®a-i,ClT,Af), (20) 
where 

Or = Si(ai ® CTi - CT 2 ® cr 2 )5i, 

A+ = 5,(0-1 ®ct 2 + ct 2 ®cri)5j, (21) 



and <5>i = & ® &, with 



Si = ST 1 



(22) 



(23) 



and 



A?=X>y) 2 > *i = V 2 & («« + #*)■ ( 24 ) 



Finally, in this case 



0W> 



iff o-j = &l; 



(25) 



span(I®I), otherwise. 



If di (or ct^ is singular, the previous computations are 
not longer valid. In this case, the commutants must be 



evaluated by direct computation and it is not possible, 
in general, to express their structure in a compact form. 

We have all the ingredients to evaluate the contribution 
related to the dissipative generators Vi in pip in every 
situation. The case in which all the A; vanish but one 
has been discussed above. The remaining cases can be 
completely described by considering the following prop- 
erties. 



(i) If A+ ^ for several indices i, 
p| {Vi,V?}' =span(I< 



(ii) If A^ ^ for several indices i, 

f| i V i> V i}' =span(I®I). 



(26) 



(27) 



(iii) If Aj = A+ 7^ and Xj = A^_ 3 ^ for a pair of 
indices («, j), then 

R^'nfe^}^ (28) 

span(I ® I, I ® (J, , (jj ® I, di ® crj) , if ct; = ct] = &j ; 



span(I ® I), 



otherwise. 



To begin with, we assume H$ — and we focus on the 
dissipative contribution to the dynamics. We denote by 
n + and rt_ the number of non-vanishing eigenvalues of 
the type A+ and A~ respectively. The relevant algebras 
in the non-trivial cases are summarized in Table HI where 
the projectors II are defined as 

n* = [7i$], n%j = 6ik8 jk , k = l,..A; 

n_ = -(i® i- n + ), n + = i®i-n_. (29) 

We notice that A — B is equivalent to n_ = 0. We 
further observe that, in case (iii), [A, B] = 0, thus it is 
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possible to choose U — U. Moreover, B = aA implies 
ctj = ctj for the index £ such that 7^ and 7^ 0. 

The sets of projectors defined in Theorem [T] are re- 
ported in the table. In the remaining cases A4 = 
span (I ® 1), part 1 of Theorem Q] applies and the sta- 
tionary state is unique. Therefore, the cases described in 
Table U are necessary conditions for the indirect manipu- 
lation of the asymptotic state of the target system T via 
the auxiliary system P. 



III. STATIONARY STATES 

We separately explore the non-trivial cases described 
in Section|TTJ Following TheoremQ]and considering Table 
HI it is possible to find the family of stationary states p^ , 
and then to extract the corresponding stationary state of 
the target subsystem, 

pF=TtApf, (30) 

by a partial trace over the degrees of freedom of the aux- 
iliary system. We consider two different choices for the 
initial state ps(0). In the spirit of the indirect control 
scheme, it can be a factor state, 



is given by 

PT = u ii (pi u S,i ~ P2 U £,2 + P3 u&) 

Pf = -u^ip^u^x- f%u& + P3V3J (34) 

PT = w ?3 (plu(i ~ P2H2 + P3H3) > 

where k = 1,2,3, are the Block vector components of 
p^?. They depend solely on the initial state of T, there- 
fore, in this case, it is not possible to manipulate p^P 
by means of an initially uncorrelated auxiliary system. 
However, if the initial state (|33[) is taken into account, a 
dependence is exhibited in terms of the Schmidt coeffi- 
cient P: 

p™ = (2P-l)u^ 3 

p? = -(2P-l)u ( 2H 3 (35) 

pf = (2P-l)uf 3 . 

Therefore, if it is possible to vary the initial degree of en- 
tanglement, different stationary states for the dynamics 
result. 



B. Case II 



Ps(0) = pt(0)® P a(0), (31) 

where pr(0) and pa(0) are arbitrary states for the two 
subsystems, that will be written using a Bloch vector 
representation as 

1 3 

Pt(0) = -(i + 5>^W), (32) 
fe=i 

and analogously for pa(0), with real coefficients p\ and 
p^. This situation refers to initially uncorrelated sys- 
tems, that will in general couple during their joint evo- 
lution. Alternatively, we consider the pure initial state 

ps(0) = |V) = Vp\ t>®| T) + Vi — p\ l)<8>\ I), 

(33) 

where Pel, and | T)> I 4) are the +1, respectively — 1 
eigenvalues of the operator er 3 . This state is entangled 
if P 7^ 0,1, and it is maximally entangled if P = i. It 
is not an arbitrary entangled state, nevertheless it can 
show the impact of initial correlations between the two 
parties on the manipulation of the stationary state of T. 

After discussing the four non-trivial cases presented in 
Table HI we recover the case with Hs 7^ 0. 



A. Case I 

In this case, A^ 7^ for some £ £ {1,2,3}. For the 
initial state (|3ip , the stationary state of the target system 



In order to fully characterize the stationary states of 
the dynamics, following Theorem [1] we need to find a 
stationary state po of maximal rank, such that 

npoy e + -i{^V e ,p o } = 0, (36) 

where £ € {1, 2, 3} satisfies Aj 7^ 0. However, it turns out 
that every solution of (f3"6"| has at least a null eigenvalue. 
In order to prove this, it is not restrictive to assume that 
only two coefficients u^ k , k = 1,2,3, are non- vanishing. 
In fact, the general case reduces to this one by a unitary 
transformation of the Pauli matrices Uk ■ Since the order 
is not relevant, we assume that 

u* a = e i/3l cos7, u\ 2 = e l/32 sin 7 , u| 3 = 0, (37) 

with real /3i, 02 and 7. The general solution of ([55]) is 

PO = i(l®I + ri(I<g>cr 3 + (7 3 <g>II)+r2Cri<X>cri (38) 

+ r 3 ((7i <g> <7 2 + 02 ® fi) + r 4 <7 2 ® (72 + r 5 cr 3 ® cr 3 ) , 

where <2j, i — 1,...,5 are real, dependent coefficients, 
such that 

ri = -i(4r 5 + l)sin(/3i-/3 2 )sin27, 

r 2 = i(4r B -l-(4rfe + l)cos27), (39) 

r 3 = -i(4r 5 + l)cos(/3i-/3 2 )sin27, 

r 4 = g ( 4r 5 - 1 + ( 4r 5 + 1) cos 2jj . 
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An explicit computation proves that Detpo = irrespec- 
tive of r 5 , and then, in full generality, there is not a 
maximal rank stationary state, in this case. Therefore, 
it is not possible to apply Theorem [TJ 



Case III 



In this case po is the solution of 

V ePo V £ + V nPo V v - \{V 2 + V 2 , po} = 



(40) 



where £ 6 {1, 2, 3}, r\ = £ + 3 are such that Ag, A J; ^ 0, 
and = V ( , V] = V n . Since 



(41) 



Vk = V A «Mf^e(I® 0- 3 + <7 3 <g)I)7£ S , 



Xr/fl^TZ^ (I ® cr 3 — cr 3 ® 1)7^ , 



the general stationary state is found to be 
/?o = -(l®I+r 1 I®a s + r-2a 6 ®l + r s a 6 <S)a i ^, (42) 

where rj, i = 1,2,3 are real, independent coefficients, 
and n = r 2 if A^ 7^ A, ( . The simplest choice leading 
to a maximal rank stationary state is = for all i. 
From it, following Theorem [TJ it is possible to build the 
complete set of stationary states of the dynamics, and 
finally to extract the stationary states of T. Although 
the algebraic structures are different, it turns out that 
the results are the same of Case I, expressed in (|34|) and 

(25). 



D. Case IV 

In this case, the maximal rank stationary state can be 
found by solving 



E(^ v *-50tfiW)=o, 



(43) 



where the 14 operators correspond to eigenvalues A^ 7^ 0, 
and k assumes two or three values in the set {1,2,3}. 
First of all, we notice that A 7^ A T is a necessary con- 
dition for the indirect control of the asymptotic state of 
T. In fact, in A = A T , all the relevant Vu operators are 
Hermitian, and then it is possible to choose the maximal 
rank stationary state as po = <B> I, leading to = 
for all k, for both correlated or uncorrelated initial states. 

For A 7^ A T , the general expression of the maximal 
rank stationary state po is rather involved, therefore we 
prefer to exhibit a concrete example proving that, in this 
case, both uncorrelated and correlated initial states allow 
indirect manipulations of the asymptotic states of T. We 
thus assume a Kossakowski matrix A of the form 



.4 



a id 
-id b 
0c 



(44) 



where a, b, c and d are real parameters satisfying the con- 
ditions 



a^O, 6^0, 0, 
ab-d 2 ^ 0, 



(45) 



expressing the complete positivity of the evolution. In 
this case, the maximal rank stationary state is found to 
be 

po = j(l®I + ri(I®(T3 + a-3®I) (46) 
+ r 2 ((Ti <E) <7i - cr 2 ® 02) + r 3 a 3 ® 03 



with 



''2 



2d 

a + b 



(b - a)d 2 



(a + b) (ab + ac + be) 

(a + b + 4c)d 2 
(a + b) (ab + ac + be) 



(47) 



The asymptotic state of T for the uncorrelated initial 
state has components 



p? = 0, p? = 0, 



(48) 



3 + 2r 2 + r 3 



(3 + £/«); 



fc=l 



for the correlated initial state we get 
p? = 0, p? = 0, 

3 + 2r 2 + r 3 V / 



(49) 



Therefore, in both cases it is possible to manipulate pjj?. 

Under the assumption Hs = 0, the only candidates 
for the realization of the asymptotic control protocol by 
means of an auxiliary system are the dynamical systems 
satisfying the conditions expressed in Case IV, with the 
additional request A 7^ A T . We now discuss the impact 
of a non- vanishing Hs ■ 

According to (fTl"]) . when adding the Hamiltonian term, 
the algebras M. computed in Section|TTJare left unchanged 
or reduced, depending on the form of Hg. In particu- 
lar, in Cases II and IV, they are not modified whenever 
[Hs, ^ + ] = 0, that is Hs is invariant under the exchange 
T <-> P (it contains only terms of the form I eg) <Ji + Ui % I 
or crj <g) crj + <jj ® Ui, with i, j — 1, 2, 3). Otherwise, M. is 
one-dimensional and part 1 of Theorem [TJ applies. 

Since in Case III the maximal rank stationary state 
previously considered is a stationary state without refer- 
ence of Hs, in Cases I and III the results obtained under a 
purely dissipative dynamics are valid in general. In Cases 
II and IV, po can be evaluated only after specifying Hs, 
therefore it is not possible to give general results. 
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IV. DISCUSSION AND CONCLUSIONS 

In this work, we have studied to what extent it is pos- 
sible to drive the asymptotic states of a system T via 
an auxiliary system P, following the indirect control ap- 
proach, when both target and auxiliary systems are two- 
level systems. We have initially considered the case of a 
purely dissipative dynamics, and then we have general- 
ized our results in the presence of an Hamiltonian term. 
We have assumed that the two systems interact sepa- 
rately with a common, homogeneous environment, that 
is invariant under spatial translation. This choice is ex- 
pressed by a particular form of the Kossakowski matrix 
C for the composite system. We have found necessary 
conditions for the indirect manipulation of the station- 
ary state of T through the initial state of P when the 
initial state is a product state, and provided a concrete 
example in which this dependence is indeed apparent. 
We have also considered the impact of an initial entan- 
glement between the two systems for control purposes. 

We have found that this kind of control can be per- 
formed only when the blocks of the Kossakowski matrix 
satisfy the condition A = B. This is not merely a mathe- 
matical request; in fact, physical system that are well de- 
scribed by a quantum dynamical semigroup in this form 
can be found in concrete experimental situationSjfor ex- 
ample in the study of the resonance fluorescence [27], HH , 
or in the analysis of the weak coupling of two atoms to an 
external quantum field [2^ |. when the distance between 
the two atoms can be neglected. 

We observe that the phenomenon described in this 
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